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Abstract 

We present a derivation of the sharp-interface limit of a generic fluctuating phase-field model 
for solidification. As a main result, we obtain a sharp-interface projection which presents noise 
terms in both the diffusion equation and in the moving boundary conditions. The presented 
procedure does not rely on the fluctuation-dissipation theorem, and can therefore be applied to 
account for both internal and external fluctuations in either variational or non-variational phase- 
field formulations. In particular, it can be used to introduce thermodynamical fluctuations in 
non-variational formulations of the phase-field model, which permit to reach better computational 
efficiency and provide more flexibility for describing some features of specific physical situations. 
This opens the possibility of performing quantitative phase-field simulations in crystal growth while 
accounting for the proper fluctuations of the system. 

PACS numbers: 81.10.Aj,81.30.Fb,05.40.-a,68.08.-p,64.70.Dv 



I. INTRODUCTION 



In recent years, phase-field models have emerged as an efficient technique to simulate 
interfacial phenomena in non-equilibrium systems [jj. This method has mainly been devel- 
oped for solidification j^, S, Q, but has also successfully been applied to other problems, 

nT n n 

such as grain boundaries |5[ , crack propagation p , viscous fingering |7J or vesicle dynamics 
^j]. The phase-field approach introduces an equation for a continuous variable <p(r, t), which 
appears as an order parameter, and takes distinct, constant values in the different phases. 
The interface is then described by the level set <ft — constant, and the transition between 
both phases takes place in a diffuse interface of thickness W. The model is completed by 
coupling the <ft equation with a diffusion field which acts as a driving force for the motion 
of the front. The behavior of the diffuse interface can then be computed by the integration 
of a set of partial differential equations for the whole system, therefore avoiding the explicit 
tracking of the interface position. This has practical advantages over using the free bound- 
ary conditions that are characteristic of a moving boundary description. Phase-field models 
are usually constructed to recover the classical moving boundary dynamics in the so called 
sharp- interface limit as W — > [9|. This limit is taken by means of a systematic asymptotic 
expansion on the interface width, and allows the model parameters to be determined in 
terms of the physical properties of the system. 

In early phase-field formulations, the model equations were derived from the variational 
minimization of a global free-energy functional for the heterogeneous system. Such varia- 
tional formulations, however, in spite of their appealing structure, presented poor compu- 
tational efficiency and did not permit to obtain truly quantitative results. For this reason, 
recently proposed phase-field models are not derived from a variational principle, but are 
specifically constructed to recover a certain moving boundary problem in the sharp-interface 
limit Besides presenting a better computational behavior, non- variational phase- field 

formulations provide for more flexibility in the description of some particular features such 
as different transport properties in the solid and liquid phases 
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ways been an important issue in the 
13, Q|. Indeed, internal or external 



On the other hand, the presence of fluctuations has a 
study of pattern-forming instabilities in crystal growth 
noises pl ay t he role of an initiation mechanism for the morphological deformations of the 
interface Thermal or solute fluctuations, for instance, must be taken into account 
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in order to study important problems such as the dynamical selection of the primary spacing 
in directional solidification [16] or the formation of secondary instabilities (sidebranches) in 
dendritic growth Fluctuations were soon introduced joto phase-field models in an ad 



hoc way as a controlled source of interfacial perturbations 



181 ] . However, phase field models 



accounting for internal thermodynamical fluctuations have not been proposed until recently, 



and in the context of variational formulations 
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22 1 . In such variational cases, the 



statistical properties of the fluctuating terms can straightforwardly be determined by using 
the fluctuation-dissipation theorem, following the lines applied by Hohenberg and Halperin 
within the context of critical dynamics j^J. In non- variational formulations, however, the 
fluctuation-dissipation relation becomes useless for this purpose because the dynamics of the 
system cannot be derived from a single free-energy functional. 

The aim of this work is to present a systematic procedure to account for the introduction 
of generic sources of noise in either variational or non-variational phase-field models. To this 
end, we will perform the sharp-interface limit of a fluctuating phase-field model for solidi- 
fication and explicitely obtain the properties of the projected noise terms that will appear 
in the moving boundary equations. This projection, which does not rely on the fluctuation- 
dissipation theorem, will be carried out by means of a hybrid asymptotic expansion which 
combines a standard sharp-interface limit with a small noise assumption for the intensities 
of the noise terms in the model. The structure of the resulting sharp-interface projection 
takes the form of a moving boundary problem, which now includes bulk and interfacial 
stochastic terms. The statistical properties of these new terms are related to those of the 
noises appearing in the starting phase-field equations. The extension of our procedure to 
thin-interface asymptotics j^] is straightforward and is not presented here for the sake of 
clarity. 

As a particular case, this analytical technique will enable a prescription for the introduc- 
tion of internal thermodynamical fluctuations in non-variational phase-field models, subject 
only to the constraint of providing the correct interface equilibrium fluctuations. This ap- 
proach will also allow for the consideration of more general noise sources of an external 
origin, such as experimental imperfections or controlled perturbations, which do not follow 
equilibrium statistics. It is worth pointing out that while the calculations will be performed 
within the framework of the symmetric solidification model, the approach can be easily 
extended to one-sided formulations jl^ . 
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This work has been organized as follows: The stochastic model equations are presented 
in Sec. HU The asymptotic stochastic procedure is developed in Sec. IHIl which has been 
divided in four different subsections: Sec. 1111 Al and 1111 Bl are dedicated to find solutions 
of the equations in the inner and outer asymptotic regions, respectively. The solvability 
conditions for the inner expansion are imposed in Sec. 1111 CI whereas in Sec. 1111 Dl we perform 
the asymptotic matching between the inner and outer stochastic fields in order to obtain the 
form of the projected equations. The projected problem is then compared in Sec. IIVI with 
the standard Lanvegin formulation for solidification allowing for the determination 

of the model parameters in the case of having internal noises of a thermodynamical origin. 
A numerical test for the validity of the approach is reported in Sec. El and Sec. |VT]is devoted 
to present some discussion and concluding remarks. 



II. MODEL EQUATIONS 

Our approach starts from a generic non-variational phase-field model, which applies for 
both the solidification of a pure substance and for the symmetric solidification of a dilute 
alloy with a constant miscibility gap ^], 

ae 2 d t (j) = e 2 V 2 (f) - /'(</>) - e\g'((f))u + ^(r, t) (2.1) 
d t u = W 2 u + ^d t h{<f>) - V • q(r, t) , (2.2) 

where a is a parameter determining the time scale of the phase-field dynamics and A accounts 
for the coupling strength between <ft and the diffusion field u. We choose g(<f)) and h((j)) to be 
odd polynomial functions of satisfying the limiting conditions g'(±l) = and h(±l) = ±1, 
and f((p) to be given by the standard double- well potential 

m = -/ - \<p 2 - (2.3) 

In the model equations Eqs. (|2.1jh ()2.2|) . u is a reduced diffusive field defined by u = 
(T — Tm)/(L/c) in the case of pure substances and by u = (c — |(c^. + c° L ))/Ac Q + \g{(f)) 
for symmetric alloys, where Tm is the melting temperature, L the latent heat per unit 
volume, c the specific heat per unit volume and Aco = c° L — c° s , being Cg, c° L the solid and 
liquid equilibrium concentrations of the alloy, respectively. The two minima <fi = ±1 of f(<p) 
in Eq. ()2.3|) correspond respectively to the solid and liquid phases of the system, so the 



interface will be represented by the transition zone between these two values. Space and 
times in Eqs. (12. lj) . ()2.2|) have been scaled out using a characteristic length I and a time 
scale v = l 2 /D, where D is the thermal or chemical diffusivity of the substance. The control 
parameter e = W/l is the scaled interface thickness, and will be the small parameter in 
which the formal expansions will be carried out. 

Fluctuations appear in the model as a non-conserved noise term 77 in the equation for 
the phase-field, and as a conserved stochastic current q in the diffusion equation. These 
fluctuations account for generic noise sources of either an internal or an external origin. We 
assume that the noises are white and Gaussian with correlations given by 



In the proposed phase field model, parameters such as a, A, and the noise amplitudes cx^, 
(j u , are intended to represent (or to be directly related to) physical parameters. On the 
contrary, the scaled interface width e has been introduced as an expansion parameter. As 
a matter of fact, Eqs. (|2.1j) . ()2.2)1 have been constructed so that the resulting dynamics (in 
the double limit of sharp interface and small noise) will be independent of e. In particular, 
the scaling factor e 3//2 of the the noise term in Eq. (|2.1|) has been introduced in order to 
make the fluctuations of the interfacial dynamics, as will be obtained below, independent of 
e. The details of this calculation and the presentation of the results are given in the next 
section. 

III. HYBRID ASYMPTOTIC EXPANSION 

In order to deal with fluctuating phase-field models, the standard asymptotic expansion, 
performed in terms of a small interface thickness, should be complemented with a small 
noise assumption. The combination of these approaches will give rise to a hybrid asymptotic 
procedure. To this end, the small noise assumption will be imposed by assuming that a^, a u 
obey some order relations with the interface thickness e. Namely, we take 




(2.4) 




(2.5) 




(3.1) 




(3.2) 
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which will permit along the expansion procedure to maintain the fluctuating terms as small 
perturbations at the desired order in a consistent way. Relations 1)3.1)1 . ()3.2)1 should not 
be understood as any explicit dependence of these parameters on e, but only as a way to 
formalize a double expansion in terms of a single vanishing parameter, the interface thickness, 
e -> 0. 

Our method closely follows the standard asymptotic procedure described in Ref. {jLf . We 
start by dividing the system into two different regions: an outer region far from the interface 
at distances much greater than e, where the phase field presents the two constant values 
= ±1 representing the solid and liquid phases at each side of the interface, and an inner 
region located around the interface up to distances of order e, where the phase field varies 
between these two values. In the limit e —>■ 0, solutions for the fields in both regions should 
match order by order in e at some intermediate distance tm, which can be taken of order 
r M ~ e 1/2 . 



A. Outer region 

In the outer region, the equations can be solved at each order by expanding the fields in 
powers of e as 

u = Uo + £Ul + £ 2 u 2 + 0(e 3 ), (3.3) 
= 0o + e0i + e 2 2 + 0(e 3 ), (3.4) 

and by expanding in Taylor series around = 0o the functions /, g appearing in the model 
equations. If we use the order relations Eqs. 1)3.1)1 . ()3.2j) . the noise terms can be assumed to 
be of orders 

r] ~ 0(e 3/2 ), (3.5) 
q ~ 0{e 2 ), (3.6) 

and the outer equations can then be obtained at each order in e. 
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1. Zero Order 



At the leading order (e°), the outer equations are given by 

/'(&) = 0, (3.7) 
d t u = V 2 « o + ^ t /i(0 o ). (3.8) 

Introducing the function /(</>) into Eq. (|3.7|) . we obtain </> = ±1, and using that h(±l) = ±1, 
Eq. (j3.8|) adopts the form 

d t u = V 2 u . (3.9) 



5. First Order 



At first order in e, we find 



r(0 o )0i = -V(0oK (3.10) 

From Eq. ()3.1U|) we determine <f>\ = by noting that the functions /, g satisfy /"(±1) 7^ 
and <?'(±1) = 0, and introducing 4>i = into Eq. (j3.11|) we get 

d t ui = VV- (3.12) 



3. Second Order 

At second order (e 2 ), and using that <pi = 0, the random current q appears in the equation 
for the outer diffusive field 

f(0o)02 = -A</(0oK (3.13) 
ftu 2 = V 2 m 2 + ^d t [h'{(f>o)(f>2) - V • q . (3.14) 

Using that </(±l) = and /"(±1) 7^ 0, equation (j3.13|) is solved by 02 = 0, and the second 
term at the right hand side of Eq. ()3.14|) can be neglected. Collecting the results obtained 
at the three first orders, the outer fields are given, up to second order in e, by 

= ±1 + 0(e 3 ), (3.15) 
d t u = V 2 m- V-q(r,t) + 0(e 3 ) . (3.16) 
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B. Inner region 



For the inner region, we write Eqs. (|2.1j) . ()2.2j) in a curvilinear coordinate system centered 
at the interface. The idea is that the solvability condition for the very existence of solutions 
of these transformed equations will provide the evolution of the coordinate system, i.e. of the 
interface, which in fact constitutes the solution we are looking for. To define this coordinate 
system by maintaining it smooth at small scales, we use an auxiliary coarse grained field 
defined as a local spatial and temporal average of the fluctuating field 0. The surface 
corresponding to the level set of this coarse grained field (<p(r,t)) = allows to define the 
3D orthogonal curvilinear coordinate system (r, si,S2), where r is a normal distance from 
the surface and Si, s 2 are the arclength distances measured along the principal curvature 
directions of the surface. Furthermore, we introduce in the inner region the scaled normal 
coordinate p = r/e and the scaled time r = t/e. We use capital letters to refer to all the 
fields when written in the inner region. After some manipulation, and keeping terms up to 
second order in e, we obtain the inner equations in the frame of the moving interface 



ae[— - vd p }<$> = d 2 <S> + end^ - e 2 p(n 2 - 2Tl)d p $ 

UT 

+ e 2 J2 - /'(*) - sXg'mU + e^Hip, s, r), 

8=1,2 



(3.17) 



vd p ]U = ±d 2 U + 1 -[k- s P (k 2 - 2U)]d p U 

(3 ' 18) 

4=1,2 

where v = v(s,t) is the local normal velocity of the interface, and we have introduced 
k = K\ + «2 and II = as the mean and Gaussian curvatures of the surface, being 

Ati(s, r), k 2 (s, r) its principal curvatures. 

The fluctuating functions H(p,s,r) = er)(r,t) and Q(p, s, r) = eq(r, t) in Eqs. (|3.17jh 
(13.18(1 stand for the renormalized noises in the inner region, and Q p corresponds to the 
normal component of the stochastic current Q. The correlations of these noise terms are 
given by 

(H(p, s, r)H(p' } s', t')> = 2a 2 5(p - p')8(s - s')5(r - r'), (3.19) 
(Q,(p,s,r)g j (p / ,s',r / )) = 2^6^ - p')8{s - s')8{r - r% (3.20) 



so that the orders in e of H(p, s, r) and Q(p, s, r) are those of and cr u respectively. Note 
that the renormalization of the noise terms is a direct consequence of the scaling of the t, r 
coordinates in the inner region. Indeed, noise terms give rise to an e factor when written in 
the inner region due to the rescaling in both time and normal distances of the delta functions 
S(p) = e6(r) and 5(r) = s5(t). 

Now we can see how the small noise assumption has been implemented in our approach. 
With the choice given by Eqs. (j3.1|) . (|3.2|) for the orders in e of the noise amplitudes, e 1 ^ 2 !! 
is 0(e 2 ) in Eq. (|3.17|) and —e~ 2 d p Q p is 0(e°) in Eq. ()3.18|) . i.e. one order higher than the 
temporal derivatives in these equations. In other words, both noise terms are first order 
perturbations for the dynamics in their respective inner equations. 

At this point, we proceed as in the outer region by expanding the inner fields and pa- 
rameters in powers of e 

U = U + eU x + e 2 U 2 + 0(e 3 ), (3.21) 
$ = $ + £$i + £ 2 $2 + O(e 3 ), (3.22) 
Ki = K i0 + SKii + 0{e 2 ) , i = l,2 (3.23) 
n = ITo + rflx + CKe 2 ), (3.24) 
v = Vo + evr + Oie 2 ), (3.25) 

and inserting the expansions into the inner equations Eqs. (J3.17j) . (J3.18)) . The inner solutions 
will be obtained by matching with the outer solutions for p — > ±oo and r — > ± , respectively. 
In the phase field equations, direct matching with the outer 0« solutions Eq. (j3.15|) provides 
the limiting boundary conditions for the terms of the inner expansion 

$o(p -> ±oo) = Tl, (3.26) 
<$>i(p -> ±oo) = 0, for i = 1,2. (3.27) 

Similarly, the matching condition for the inner diffusion field requires that, at leading order, 
the gradients of Uo vanish 

lim d p U = 0. (3.28) 

p^±oo 

At higher orders, the matching conditions for the diffusive field present some additional 
difficulties due to the apparition of random terms, and will be discussed in detail in Sec. IIIIDI 
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1. Zero order 

At leading order (e° for the $ equation, e~ 2 for the U equation), the inner equations are 
given by 

d 2 p % - /'($„) = 0, (3.29) 
d 2 p U = 0. (3.30) 

Inserting the double- well potential Eq. ()2.3|) into Eq. (|3.29|) . we obtain the standard kink 
solution for the phase-field at zero-order 



$o(p) = - 




(3.31) 



which satisfies the matching condition Eq. (J3.26JI for p — > ±oo. Integrating Eq. ()3.30|) twice 
over p, we have 

U (p,8,r) = A(s,r) + B(s,r)p, (3.32) 

where A and B are integration constants. Imposing the matching condition Eq. ()3.28j) . we 
determine B(s, r) = and obtain a p-independent solution for the diffusion field at zero 
order 

U (s,t)=A(s,t). (3.33) 



2. First Order 

Using the solutions obtained at zero-order, the first-order inner equations (e 1 for the $ 
equation, e -1 for the U equation) read 



= -(voa + ^dpQo + Ag'i^Uo, (3.34) 
= ^ + ^($ ), (3.35) 



where <1 is the self-adjoint operator « . « - /*(*»,) and we have used that d^/dr = 
from Eq. ()3.31|) . As described by Almgren [11| . an expression for $! can be obtained from 
Eq. ()3.34|1 by inverting the operator Q, leading to 

$x = n-'l-ivoa + k )9 p <1>o + \g\%)U ]. (3.36) 
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Since Q is an even operator and d p &o, g'($o) are even functions of p, $1 is an even function 
of p. Integrating Eq. ()3.35|) twice over p, we get 

U 1 =D( S ,r) + C(s,T)p + ^ [' dfth{*o) + (3.37) 

J 

where D and C are integration constants and we have used that 8 p Uq = (cf. Eq. ()3.33j) ). 



3. Second order 

The stochastic terms appear in the inner equations at second order (e 2 for the $ equation, 
e° for the U equation), which are given by 

fi$ 2 = - (avt + /ti)<9„$o - (av + K, )d p §i + a-^ 

dr 

+ lf"(^l + p(4 -2U )d p % ( 3 - 38 ) 
+ A(/($ )tA + \g"($o)$iU - e- 3/2 H, 

d 2 p U 2 = - (v + K W + ^ + ^d p h($ ) + ?0p[/i'(*o)*i] 
r or I 1 

-E^o-^(*o)fi + >Q, <M9) 

i=l,2 

The first equation Eq. ()3.38|) will be used in the next section when imposing the second 
order solvability condition of the problem. Integrating Eq. (j3.39J) twice over p, we find 

U 2 = F(s, r) + E(s, r)p - (v + k ) [ dp'U x 

Jo 

+ J' dp' jT dp"^± - \d 2 s U,p 2 + | jT dp'h(%) (3-40) 

+^ [ p dp'ti(%)$ 1+ \ r dp'Q p { P ,s,T), 

A Jo £ Jo 

where F and E are again p-independent integration constants. 



C. Solvability conditions 

We impose now the solvability conditions for the inner problem, which at first and second 
orders are respectively given by 

/oo 
(d p $o)n$ j dp = 0, for j = 1, 2. (3.41) 
-oo 



Inserting Eq. (|3.34|) into the first order solvability condition, we get 



-(oeuo + Ko)h - XhU = 0, (3.42) 
which allows to determine Uq as 

U (s, r) = ~^v - t^k , (3.43) 
where I\ and li are new integral constants given by 



h = / dp(d p $ ) 2 , (3.44) 

J -co 

/CO 
dpg'{%){d p %). (3.45) 
-co 

Imposing the second order solvability condition Eq. (j3.41|) . and taking into account the 
parity of the potentials f,g,h and of the inner solutions $0, ^i, we determine an expression 
for the constant D in Eq. ()3.37)1 

D[s , T) = _ (a „, + Kl) Jl + „„ A + IL + 2k _ e -vM (,«) 
where ^3, I4 and 1$ are defined by 

/CO t-p 
dp («9 P $ ) g'($ ) / dp'h($ ), (3.47) 
■00 Jo 

ATT poo 

h = J dp {d p %)g'{%)p\ (3.48) 
dp(d p %)^, (3.49) 

and Z is a stochastic term given by 

/CO 
dp(d p Q )H(p,s,T), (3.50) 
-co 

whose statistical properties are given by 

(Z(s, r)Z(s', r')> = 2/ l( 7^(s - s')5(r - r'). (3.51) 

D. Matching of fluctuating fields 

At this point, we continue by imposing the remaining asymptotic matching conditions of 
the problem. However, the matching of the diffusion field presents some subtleties due to its 
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fluctuating character. The main problem is that, at second order in e, the U field fluctuates 
in the normal direction (cf. Eq. ([3.40)) ) . and hence cannot be written as a simple asymptotic 
expansion for p — > ±00, preventing the matching with the outer field. This difficulty can be 
overcome by introducing an auxiliary matching function defined in both regions as 

X(r,s,t) =u(r,s,t)- / dr'q r {r', s', *), (3.52) 

Jo 

X(p,s,r) = U(p,s,r)- [ P dp'Q P (p',s } r). (3.53) 

Jo 

In view of Eq. (|3.40|) . it is easy to see that the inner auxiliary function X introduced 
in Eq. (J3.53j) is smooth up to order e 2 in the matching region tm- Explicitely, if X is 
asymptotically expanded for p — > ±00 as 

X ~T + S p + Rp 2 + 0(p 3 ), (3.54) 

and the outer matching function \ is expanded in Taylor around r = O 1 * 1 by 

X ~ X^) + 9 rX (0 ± ) • r + i^ 2 x(0 ± ) • r 2 + 0(r 3 ), (3.55) 

the inner and outer terms can be matched at tm ~ £ 1//2 in the limit e — > to obtain the 
matching relations 

T = X (0 ± ), (3.56) 
S = edrxiO^), (3.57) 
R = y^ 2 x(0 ± ). (3.58) 

The last step is to expand the previous Equations (|3.56|) . (|3.57|) and (|3.58j) in powers of e to 
complete the matching at each order in e. 



1. First Order 



At first order in e, the inner field U\ given by Eq. ()3.37j) can be asymptotically expanded 
for p — > ±00 as 



TT ldU 2 , 



2 dr 

where D is given by Eq. (|3.46|) and 



P + D+ V fj±, 



±00 



dp[h(%(p))±l], 



(3.59) 



(3.60) 
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where we have used that h satisfies h(±l) = ±1 and the far field condition Eq. (|3.26|) . Since 
h(<&o) is an odd function of p, we have 

J+ = Jf = J x . (3.61) 



2. Second Order 

Similarly, the second order inner solution for the diffusive field Eq. ()3.40|) can be expanded 
asymptotically for p — > ±oo as 

cp dTL CP CP 1 dUi 
U 2 „ F + E p - (v + «o) J d P'^ + J o dP I dp"^ 

*u v X 

+ + y^i T yP + ^ y dp'Q p (p,s,t), 



where 



p 



J±= / dp'hX^o)®!, (3.63) 

and we have used the far field conditions Eqs. ()3.26|) and ()3.27|) . Inserting the expressions 
Eqs. (|3.43|) . (|3.59|) and (|3.62j) into the right hand side of Eq. (|3.53|) . we can determine the 
parameters R, S and T in the far field expansion of the matching function X (cf. Eq. (|3.54j) ) 
and perform the matching with the outer function x (cf. Eq. (|3.55p ). 

Imposing the third matching condition Eq. ()3.58|) at first order in e, we determine that 

^ = 0, (3.64) 

GST 

which, using the relation Eq. (|3.36|) . brings to 

— = 0, (3.65) 
and therefore the integral constants I4 and I5 defined in Eqs. (|3.48|) . ()3.49|) vanish 



J 4 = J 5 = 0. (3.66) 

From the two first orders of Eq. (J3.56)) , we get an expression for the outer diffusive field at 
the interface valid up to first order 

u(0±) = -±(av + «) + f^(| + J,) _ Z ±L<1 + o{e% (3.67) 
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where z(s,t) = Z(s,r)e 1 ^ 2 is a stochastic term whose statistical properties can be deter- 
mined from Eq. (j3.51j) and are given by 

(z(s, t)z(s r , if)) = 2hal8(s - s')5(t - if). (3.68) 

Note that the projected interfacial noise term has neither in Eq. (jH.fi 7j) nor in Eq. (j3.68j) any 
explicit dependence in e, which is a direct consequence of the e 3//2 factor introduced in the 
noise term of Eq. (j2.1j) . Indeed, this is the reason why such factor was introduced in the 
formulation of the model. 

The calculation is completed by imposing the matching condition Eq. (j3.57j) up to second 
order, which can be written as 

vo + ev! = d rX \+ + 0(e 2 ). (3.69) 

Inserting Eq. (j3.52j) into Eq. (j3.69jl . we get a heat /mass conservation equation valid up to 
first order in e, 

v = VQ + £Vl = [ 9rU ]- _ [ Qr ]- + (e 2 ), (3.70) 

where q r accounts for a normal stochastic current across the interface. This term, although 
being of order e 2 , has not been neglected in Eq. (j3.70j) in order to not break mass conservation 
in the stochastic diffusion equation Eq. (j3.16jl . which is valid up to second order. 

This last equation completes the sharp-interface projection of the stochastic phase-field 
model of Eqs. (j2.1j) - (j2.5j) . This projection constitutes the main result of this paper, and is 
given by the diffusion equation Eq. (j3.16j) . with the noise of Eq. (j2.5jl . supplemented with 
two moving boundary conditions at the interface: the conservation condition Eq. (j3.70jh 
and the Gibbs-Thomson Eq. (j3.67j) . where a new projected interfacial noise appears with 
correlation given by Eq. (j3.fi8j) . Note that the projected boundary conditions at the interface 
Eqs. (j3.7()j) . (j3.67j) are obtained at the order immediately lower than the order at which the 
asymptotic expansion is performed. 

While the general lines of the calculation follow the standard sharp-interface asymptotics, 
we have included the fluctuation terms during all the procedure, which have been projected 
in the weak noise limit. This calculation is thus similar to the front dynamics projection 



performed in Ref. |24(. Indeed, the projected interfacial noise appearing in Eq. (j3.67jl is the 



analogous counterpart of the noise term of the projected eikonal front equation of Ref. 
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IV. INTERNAL FLUCTUATIONS IN A GENERIC PHASE FIELD MODEL 



Thus far, the noises considered in this work are intended to account for both external 
and internal sources of fluctuations. Nevertheless, it is worth pointing out that the resulting 
stochastic sharp-interface equations are similar to those postulated in the Langevin formu- 



lation of solidification due to Karma 



14j (see also Ref. 3])' which was constructed to 



follow equilibrium statistics. This offers the possibility of using the results above to provide 
generic (not necessarily variational) phase-field models with the correct equilibrium fluctu- 
ations. To illustrate this, let us consider the Langevin sharp-interface equations ^| 

d t usi = VVi-V-q SI (r,i), (4.1) 
vsi = [d r u S i]+ ~ [gf (4.2) 
u SI (0) = -d K- pv + 9(r,t) , (4.3) 

where q s/ and 9 are fluctuating terms with correlations given by 

( 9 f (r,t)gf (r',t')> = - r')5(t - t'), (4.4) 

(0(*,tWJ)) = 2KB Xf^ (s-s')6(t-t'). (4.5) 

The Gibbs-Thompson equation Eq. (J4.3j) can be compared with Eq. (J3.67j) and the diffusion 
equation Eq. (|4.1|) with Eq. (|3.16j) . This comparison enables the determination of the phase- 
field parameters in terms of physical and substance parameters, which are given by the 
equations 

A = (4.6) 

<* = ■§-, (4-7) 
2 _ K B T 2 M c 

h K B TlcP 

+ d\m d " 1 ] 

In the last relations, the two first equations Eqs. (|4.6|) . (|4.7j) are the usual expressions deter- 
mined by the standard asymptotic procedure, whereas Eq. (J4.8|) reflects the identification 
between the conserved stochastic currents of both phase field model and sharp interface pro- 
jection. In this sense, a major result of our approach has been the derivation of an expression 
for the noise strength of the phase-field, that is Eq. ()4.9|) . from the above calculations. 
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With this election of the model parameters, the phase-field simulations will present the 
correct equilibrium statistics in the limit of small interface thickness e —>■ 0. Therefore, 
the non- variational phase- field formulation of Eqs. (|2.1jl . (|2.2j) can be used to quantitatively 
account for thermodynamical fluctuations in solidification processes. 

V. TEST OF THE APPROACH 

In order to test the validity of our approach, we have performed 2D phase-field simulations 
to obtain the power spectrum of the interfacial fluctuations of a solid-liquid stationary flat 
interface. Introducing the Fourier transform of the interface position £(r, t) as £fc(t) = 
f dk £(r, t) e~* kr , the power spectrum of a stationary planar front in scaled variables is 
given by 

S(k) = (te. k ) = J T^itkb) = (s- 1 ) 

where 7 = l d L 2 do/TMC is the scaled interfacial surface energy. In the simulations, space and 
times have been scaled using length and time scales of I = 10~ 8 m and v = 9 x 10 -10 s, 
respectively. The functions h, g have been chosen to be h(<p) = <j) an d g'{4>) = (1 — <fi 2 ) 2 so 
that the model does not have a variational structure. The substance parameters used in 

3 values of the pure SCN in the 3D case, and are given by 
and a\ = 0.001432. For this choice, and using Eqs. (|4.6|) - 
()4.9jl . the phase-field parameters take the values A = 3.13, a = 10.76 and = 0.05158. The 
interface thickness has been taken to be e = 0.3. 

The simulations have been implemented with a finite differences scheme on a 50 x 512 
lattice with Ax = Ay = 0.2 and At = 0.005. We have used the initial conditions 
(j)(x,y,0) = — tanh(x/£\/2), u(x, y, 0) = 0. Non-flux and periodic boundary conditions 
have been imposed in the x and y directions respectively. 

The numerical implementation of the stochastic terms has been carried out by generating 
Gaussian-distributed random numbers at each of the lattice sites. The correlations of these 
numbers can be determined by discretizing the time and spatial delta functions in Eqs. (|2.4|) . 
(12. 5|) by substituting 8(x — x') — > Sw/Ax and 8(t — t') — > 5 nn >/At. The divergence of the 
stochastic current in Eq. (|2.2|) has been discretized by using a forward differences scheme 
V • q(r, t)|ij = (q x (i + 1, j) - q x (i,j))/Ax + (q y (i,j + 1) - q y (i,j))/Ay. 

The power spectrum statistics has been obtained as a time average among the last 3 x 10 6 



the simulations correspond to t 



d Q = 0.2817, 8 = 3.0331 |2d 27 
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FIG. 1: Comparison between the theoretical power spectrum of the stationary interface and the 
results from the phase-field simulations. 

time steps in a long-term simulation of 3.5 x 10 6 steps, and is represented by a dashed line 
in Fig. n The solid line in Fig. ^ depicts the theoretical prediction given by Eq. (|5.1|) 
and, as it can be seen, an excellent agreement is found between theoretical and numerical 
results. The vertical dashed line in the figure represents the wavelength associated with the 
effective thickness of the interface, and determines the expected breakdown of the phase-field 
description. 



VI. DISCUSSION AND CONCLUSIONS 



To summarize, we have obtained an asymptotic projection of the fluctuating phase-field 
equations (|2.1|) . (j!2.2j) to a sharp interface description. This has been worked out by means of 
a hybrid asymptotic procedure, combining sharp interface and small noise limits. As a result, 
the projected equations adopt the form of a moving boundary problem with a conserved 
stochastic force q in the equations for the diffusion field Eqs. ()3.16|) and ()3.70|) . and an 
interfacial noise z(s, t) in the Gibbs-Thompson condition Eg. ()3.67|) . Other authors have 




previously introduced fluctuations in phase- field models [20|, |2l|, |22| , but their approaches 
applied only for the case of variational formulations and were restricted to noises from a 
thermodynamical origin. 

In this context, it has been claimed (^(J that the presence of a non-conserved noise such 
as r] in the equation for the phase field is not relevant for the dynamics of the phase field 
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model, and thus could be omitted in simulations. In order to check the importance of 
the non-conserved phase field noise, we have carried out a numerical test with the same 
parameters reported in section but taking = 0. In this case, the power spectrum 
is plotted as a dotted line in Fig. ^ The clear disagreement with both the theoretical 
prediction and the simulations of the complete model indicates that the phase field noise is 
indeed necessary in order to obtain quantitative results. This can be explained by noting that 
Eq. ()3.67|) establishes a direct relation between the non-conserved phase-field noise rj and 
the interfacial fluctuations appearing in the Gibbs- Thompson equation, usually associated 
to kinetic attachment effects |l3j]. Thus, the small significance of the 77 noise reported in 
Ret. Q is probably due to tbe fact that tbe stationary power speetrum was calculated in 
the limit of vanishing kinetics (3 = 0. Therefore, we conclude that, in the presence of kinetic 
effects, the phase field noise is relevant for a quantitative description of the solidification 
process. 

It is interesting to discuss the scaling of the noise terms as proposed on the one hand in 
Eqs. (JHHJ), ()3.2|) and on the other hand in the e 3 ^ 2 factor explicitly appearing in the equation 
for the phase field, Eq. (|2.1j) . As it has already been commented, the assumption of the order 
relations Eqs. (j3.1|) . (|3.2|) has permitted to manage a double expansion (sharp interface and 
small noise) by formally using a single small parameter. The specific powers of e appearing 
in these relations have been chosen for maintaining fluctuations as small perturbations for 
the dynamics of both inner and outer equations. On the contrary, the multiplicative factor 
of the noise term of the equation for the phase field, Eq. (|2.1)l . has a different motivation. It 
is well known that there are problems in the formulation of stochastic field equations when 
the noise terms are delta-correlated in space. I such cases, some kind of regularization is 



required. In Ref. 



24j , for instance, this regularization was provided by the correlation length 



of the noise, in such a way that the results did depend on that parameter. In the present 
case, the regularization is provided by the interface width e. The projection of the bulk 
noise into the interface gives a fluctuation term that in principle should diverge as e goes to 
zero. The e 3 / 2 factor of the noise term in Eq. (|2.1jl exactly cancels out this divergence, and 
has been introduced in the formulation of the model precisely to make results independent 
of e, specifically regarding the new interfacial noise term z in Eqs. (|3.67|) . (|3.68|) . 

In conclusion, we have proposed an asymptotic procedure to obtain the sharp-interface 
projection of a generic (not necessarily variational) phase-field model with fluctuations. We 



19 



have tested the validity of our approach by comparing the phase-field results with the the- 
oretical prediction for the interfacial fluctuations in a simple solidification problem. This 
procedure can be useful both in situations where the fluctuation-dissipation theorem does 
not hold, such as in presence of external fluctuations, and when the phase-field model has 
a non-variational nature. An important example of this latter case corresponds to efficient 
phase-field models for the solidification of alloys 0, 12 1, which with our method can incor- 
porate internal fluctuations. The availability of these models for quantitative simulations 
in the fluctuating case appears as a promising step towards the study of complex situations 
such as the apparition of dendritic sidebranching or the wavelen gth selection during initial 



redistribution transients in the directional solidification of alloys [2i 
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